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ABSTRACT 

We present a comprehensive study of the spectrum of the narrow-line Seyfert 1 galaxy 
RE J1034-I-396, summarizing the information obtained from the optical to X-rays with obser- 
vations from the William Herschel 4.2m Telescope (WHT), the Hubble Space Telescope (HST), 
the Extreme UltraViolet Explorer {EUVE), ROSAT, ASCA and Beppo-SAX. The Beppo-SAX 
spectra reveal a soft component which is well-represented by two blackbodies with A:Toff=60 eV 
and 160 eV, mimicking that expected from a hot, optically-thick accretion disc around a low-mass 
black hole. This is borne out by our modeling of the optical to X-ray nuclear continuum, which 
constrains the physical parameters of a NLSl for the first time. The models demonstrate that 
RE J1034-I-396 is likely to be a system with a nearly edge-on accretion disk (60 to 75° from 
the disk axis), accreting at nearly Eddington rates (0.3 to 0.7 LEdd) onto a low mass black hole 
(Mfe/j ^2 to 10x10^ M0). This is consistent with the hypothesis that NLSls are Seyfert-scale 
analogies of Galactic Black Hole Candidates. The unusually high temperature of the big blue 
bump reveals a flat power-law like continuum in the optical/UV which is consistent with an 
extrapolation to the hard X-ray power-law, and which we speculate may be similar to the con- 
tinuum component observed in BL Lac objects in their quiescent periods. From the Beppo-SAX 
and ASCA data, we find that the slope of the hard X-ray power-law depends very much on the 
form of the soft component which is assumed. For our best-fitting models, it lies somewhere 
between a=0.7 and 1.3 and thus may not be significantly softer than AGN in general. 

Subject headings: accretion, accretion disks - galaxies: active - galaxies: nuclei - galaxies: Seyfert - 
X-rays: galaxies - galaxies: individual (RE J1034-I-396) 



1. Introduction 

The narrow-line Seyfert 1 (NLSl) typo of 
AGN has proven to be a valuable and a fasci- 
nating resource for the study of the class as a 
whole. They have very strong soft X-ray excesses 
(Puchnarewicz et al. 1992; Boiler, Brandt & Fink 
1996), arc often highly variable (Brandt, Fabian 
& Pounds 1996; Boiler et al. 1997) and, in some 
cases, the optical/UV big blue bump (BBB) com- 
ponent is so hot that it is shifted into the UV/EUV 
regime, leaving a bare, power-law-like contimuim 
component in the optical (Puchnarewicz et al. 
1995). 

One of the earliest hypotheses put forward to 
explain the strong and variable ultra-soft X-ray 
excesses, was a high mass accretion rate onto 
a relatively low-mass black hole (Pounds, Done 
& Osborne 1995). This was originally proposed 
for RE J1034-I-396, by analogy with the proper- 
ties of Galactic Black Hole Candidates (GBHCs). 
An alternative model was that such systems were 
geometrically-thick accretion disks (ADs) viewed 
face-on, lying co-planar with a flattened broad line 
region [BLR; Puchnarewicz et al. (1992)]. 

In our initial study of RE J1034-I-396 (Puchnarewicz 
et al. 1995), we showed that the overall IR to X- 
ray spectrum compared well with the combination 
of a geometrically-thin, optically-thick AD and an 
underlying power-law with a slope, q;=1.3 (a is 
defined throughout such that oc z/~"). The 
Kerr black hole mass was 7x10^ M©, the accre- 
tion rate was 0.073 Mq yr~^ and the disk was 
viewed at an angle of 60° from its axis. However, 
this was for illustrative purposes only and no con- 
straints were placed on these parameters. Thus, 
while this might have suggested the presence of a 
low-mass black hole for RE J1034-I-396, consistent 
with the GBHC analogy, solutions at higher black 
hole mass could not be ruled out. 

With the launch of Beppo-SAX and the inclu- 
sion of RE J1034-I-396 in the NLSl core program, 
came the opportunity to place meaningful con- 
straints on black hole mass (Mj/j), accretion rate 
(M) and inclination. RE J1034-I-396 was also due 
to be observed with HST within a few months of 
Beppo-SAX, and an optical spectrum was also 
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scheduled in WHT service time within a few 
weeks. Thus a quasi-simultaneous spectrum, pro- 
viding the most complete coverage possible from 
the optical to mid X-rays was available. Simul- 
taneity is important for a NLSl due to the vari- 
able nature of the class, although RE J1034-I-396, 
contrary to its counterparts, has shown remark- 
able stability. Thus these data have provided 
the best opportunity yet to fit AD models and 
provide constraints on the defining parameters of 
an NLSl. Using a combination of a power- law 
and an optically-thick, geometrically-thin Kerr 
AD model, we present fits to the observed spec- 
trum of RE J1034-h396, measuring M^h and M 
and constraining the inclination for the first time, 
and discuss the results in the context of the GBHC 
model. 

The improved UV and optical data presented 
here have also allowed us to re-visit the issue of 
what is producing the optical/UV continuum in 
this object. Our previous work had suggested a 
power-law-like continuum rising towards the red, 
which was not due to the host galaxy or the BBB. 
We have separated out the galaxy component from 
the spatially-resolved optical spectrum obtained 
and, when combined with the HST data, examined 
the pure nuclear component with a greater degree 
of clarity. The results and their implications for 
the production of the optical continuum in AGN 
are also discussed in this paper. 

2. Data 

RE J 1034+396 has been observed on several 
occasions at different wavelengths from the ra- 
dio to X-rays (Puchnarewicz et al. 1995; Pounds, 
Done & Osborne 1995; Puchnarewicz, Mason & 
Siemiginowska 1998; Breeveld & Puchnarewicz 
1998). In this paper, we present data taken at 
optical, EUV (EUVE), soft X-ray {ROSAT HRI 
and SAX LEGS) and medium X-ray {SAX LEGS 
and MEGS) energies which were taken quasi- 
simultaneously with the ffST-FOS spectrum. We 
also analyze previously unreported ROSAT HRI 
(taken 1994 November 20) and ASCA data (taken 

1995 May 18) to check for variability in this object. 

2.1. Optical spectra 

RE J1034-I-396 was observed on the night of 

1996 March 24 using the twin-armed ISIS spec- 
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trograph on the William Hcrschel Telescope on La 
Palma as part of the service observing program 
(observers D. Pollacco and D. King). The input 
starlight was divided between the red and blue 
arms of the spectrograph by means of a dichroic 
filter whose crossover wavelength is 5400A. The 
R316R grating was used in the red arm of the 
spectrograph, centered on a wavelength of 6677 A. 
The blue arm utilized the R300B grating centered 
on 4519 A. The red spectrum spanned the wave- 
length range 6OOOA to 7700A, while the blue spec- 
trum extended from 3800A to 5300A. The data 
were recorded on TEK 1024 square CCDs. Two 
exposures, each of 600s duration, were taken of 
RE J 1034+396 and the resulting images combined 
during analysis (eliminating cosmic ray events) to 
give a total effective exposure time of 1200s. The 
sky was clear throughout. Standard exposures 
of wavelength calibration arcs and flux calibra- 
tion stars were made immediately following the 
RE J 1034+396 observations using the same instru- 
mental set-up. 

The integrated spectrum of RE J1034+396 is 
clearly contaminated by light from the host galaxy, 
particularly in the red. The separate contributions 
of the active nucleus and the underlying galaxy 
are clearly evident in the cross-dispersion spatial 
profile. To remove the worst effects of the galaxy 
contamination on the nuclear spectrum, the 2-D 
image of the spectrum was divided into 25 A bins 
and the mean spatial profile in the cross-dispersion 
direction formed for each bin. Two Gaussian pro- 
files superimposed on a flat background were then 
fitted simultaneously to the spatial data in each 
bin, to model the nuclear and extended compo- 
nents respectively. Gaussian functions were found 
to be satisfactory and robust representations of 
the spatial data for the purposes of this decom- 
position. The flux in the two Gaussian compo- 
nents as a function of wavelength thus represents 
the spectrum of the nuclear and extended emission 
respectively. 

2.2. HST observations 

RE J1034+396 was observed by HST on 1997 
January 31 using three gratings (G130L, G190L 
and G270L) covering the range 1100 A to 3300 A. 
Full details of the observations and data reduc- 
tion were presented in Puchnarewicz, Mason & 
Siemiginowska (1998). 



The HST spectrum was combined with the de- 
convolved optical nuclear spectrum and the re- 
sulting optical/UV continuum was fitted with a 
power-law model. The best fit was obtained for 
a=0.9, having removed obvious emission and ab- 
sorption features before fitting. Errors are difiicult 
to determine precisely, because the relative nor- 
malization of the two spectra is not known (the 
data were not taken simultaneously and difli'erent 
slit widths and positions were used) and because 
of the difficulty in securely removing absorption 
and emission features before fitting the continuum. 
Assuming, therefore, a conservative error of 20% 
on the optical and UV fiuxes, gave a greater than 
99.9% probability that the model was a good fit 
to the data. 

2.3. EUVE data 

RE J1034+396 was observed with the Deep 

Survey Spectrometer (DS/S) on board EUVE 
from 14 April 1997 UT 05:35:04 to 20 April 1997 
UT 01:03:07 for a total of approximately 100 
ks. The observation was part of a simultaneous 
EUVE/Beppo-SAX campaign. 

The DS/S (Welsh et al. 1990) is equipped with a 
broad band imaging detector (covering the 66 — 178 
A or 0.07 - 0.18 keV band in the Lexan/B fil- 
ter) and three spectrometers (Hettrick & Bowyer 
1983; Abbot et al. 1997) covering the "short" (SW: 
70 - 190 A, 0.06 - 0.18 keV), "medium" (MW: 
140 - 380 A) and "long" (LW: 280 - 760A) EUV 
wavelengths. This configuration allows simultane- 
ous imaging and spectroscopy with a spatial reso- 
lution of ~ 1 arc min and a spectral resolution of 
A/AA ~ 200 at the short wavelengths. 

RE J1034+396 was observed 0.2° off-axis in or- 
der to avoid the DS dead-spot and to extend the 
spectrum toward shorter wavelengths (Amin ~ 66 
A) with respect to a regular on-axis observations 
(Amin ~ 72 A). [The DS dead-spot is a small re- 
gion of reduced gain and detector quantum effi- 
ciency near the center of the DS detector, caused 
by the observation of the very bright EUV source 
HZ 43 (Sirk 1993).] 

2.3.1. Light- curve 

We extracted the light-curves from the DS time- 
ordered event list using the EUVE Guest Observer 
Center software (IRAF/EUV package) and other 
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IRAF timing tasks adapted for EUVE data. Af- 
ter correcting the data for instrumental dead-time, 
telemetry saturation, vignetting and eliminating 
intervals of high particle background, the total ef- 
fective exposure is 97061 s. 

We counted the source photons in a circle of 1.1' 
radius, and we estimated the background in a con- 
centric annulus with inner and outer radii of 1.5' 
and 3.8' respectively. The extraction region in- 
cludes more than 98% of the DS point spread func- 
tion (Sirk et al. 1997). The effective area of the DS 
instrument (25 cm^ at A = 85 A), is more than 10 
times larger than the spectrometer effective area 
at this wavelength and a good detection of the 
source (average signal-to-noise ratio (SNR) ^ 4.9) 
was obtained during each orbit. Figure 2 shows 
the light curve binned over one average EUVE or- 
bit (about 5544 s). Since the EUVE instruments 
are shut down during satellite daytime (around 
2/3 of the total orbital time for EUVE 550 km 
circular orbit with an inclination of 28°) and pas- 
sages through regions of high particle background 
(such as the South Atlantic Anomaly) the aver- 
age effective exposure per orbit is approximately 
1000 s, after discarding the bins with effective ex- 
posure less than 200 s. 

The average count rate over the entire observa- 
tions is 0.038 ± 0.001 counts s~^. As a compari- 
son RE J 1034-1-396 was observed twice before by 
EUVE: during the all-sky survey in 1992-1993 it 
had a count rate of 0.023 ± 0.007 in ~ 1600 s of 
observation (Fruscione 1996) and again serendip- 
itously on 27 December 1993 for an effective ex- 
posure time of ~ 5000 s; the count rate then was 
0.029 ±0.003 (Fruscione, private communication). 
During the observation described in this paper 
(April 1997), RE J1034+396 was in a relatively 
bright state, ^60% brighter than in the all-sky 
survey and ~30% brighter than in December 1993. 

As immediately visible from the light-curve, the 

EUV flux from RE J1034+396 remained effec- 
tively constant for the 5 days of observation and 
in fact a test for a constant source on the entire 
data set (exchiding intervals with effective expo- 
sure below 200 s) gives a value of 51 for 75 degrees 
of freedom, corresponding to a 1% probability 
that the source varied during the six days of this 
observation. 



2.4. ROSAT HRI images 

Since the re-discovery of the active nucleus in 
RE J1034-h396 by the ROSAT Wide Field Cam- 
era in 1990, several observations have been made 
in the soft and medium X-rays and a log of the 
available data to date is given in Table 1. 

A measurement of the soft X-ray flux from 
RE J 1034+396 was made with the HRI on ROSAT 
in 1996 November and yielded a count rate of 
0.61±0.02 count s~^ with no evidence for variabil- 
ity for the duration of the observation (~30 ks). 
The full width at half maximum (FWHM) of the 
profile of RE J1034-h396 was 7 arc sec. The point 
spread function (PSF) of an unresolved on-axis 
source observed by the HRI has a FWHM between 
~5.1 and 7.4 arc sec, thus the soft X-ray emission 
of RE J1034-h396 in the ROSAT HRI was consis- 
tent with a point source. 

A previous HRI observation of RE J1034-h396 
was made in 1994 November 20. We measure a 
count rate of 0.49±0.02 from these data, thus the 
soft X-ray flux had increased by ~20% in the two 
years since the 1994 pointing. This is consistent 
with EUVE observations which showed an increase 
between 1992/1993 and 1997 (see Section 2.3.1). 

2.5. ASCA observations 

RE J1034+396 has been observed on two oc- 
casions by ASCA, in 1994 November (coincident 
with the first ROSAT HRI observation) and six 
months later in 1995 May. The first dataset was 
analyzed by Pounds, Done & Osborne (1995) but 
the second is unpublished to date. We present the 
1995 data here and re-analyze the 1994 spectra 
to ensure a consistent approach towards the data 
reduction when comparing observations. 

Both sets of SIS observations were taken in 1- 
CCD mode. The GIS source spectra were ex- 
tracted using a circle of radius 4 arc min. The GIS 
background regions were selected from a position 
opposite that of RE J1034-I-396, i.e. at the same 
radial distance, and with a circle of radius 6.7 arc 
min. The SIS spectra were extracted using a circle 
of radius 3 arc min and a nearby circular region of 
the SIS chips was used for the background. The 
background light-curves were used to check for and 
discard any periods of anomalous data before the 
spectra were extracted. 

The background subtracted light-curves of the 
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target were examined for signs of significant vari- 
ability. None was seen in the 72020000 data [con- 
sistent with the results of Pounds, Done & Os- 
borne (1995)]. In the 72020010 observation, there 
is a suggestion of variability at the 20-30% level, 
although a test indicates that this cannot be 
distinguished from flux at a constant level. The 
count rate of RE J1034-h396 in 1995 May was con- 
sistent with that six months earlier in 1994 Novem- 
ber at the 2(7 level. 

The reduced spectra were fitted with a combi- 
nation of power-law and thermal models, using the 
XSPEC spectral fitting package (Version 11.00) and 
the results are given in Table 1. The fitting ranges 
were 1.0 to 10 keV for the GIS and 0.7 to 10 keV 
for the two SIS. In all the fits, a Galactic column 
density of 1.5x10^° cm~^, interpolated from the 
maps of Stark et al. (1992), was assumed. For the 
earlier 72020000 observation, the lowest x^s were 
obtained for the two power-law and two blackbody 
plus power-law fits (full details are given in the 
table). For all fits, the best fit intrinsic column 
density converged to 0, indicating little if any ab- 
sorption local to the AGN. The slopes of the hard 
X-ray power-law for the two best-fitting models 
were a ~0.7-0.8 which is typical of AGN in gen- 
eral (Mushotzky 1984; Comastri et al. 1992) and 
not significantly softer as suggested by Pounds, 
Done & Osborne (1995). 

The second dataset, 72020010, was observed 
for only one-third as long as the earlier one, thus 
the fits are more poorly constrained. Only a sin- 
gle power-law model gave an acceptable fit to the 
data if the intrinsic column was allowed to be 
free: for the other three models tried, the intrin- 
sic column converged to a few times lO'^^ cm^^. 
Such high column densities were in turn offset 
by un-physically strong un-absorbed components, 
suggesting some flattening of the spectrum to- 
wards lower energies. A broken power-law model 
with no intrinsic column did give an improved fit 
over a single power-law. The single power-law 
fit could not be improved by the addition of a 
bremsstrahlung or blackbody component however. 

The 72020010 data were well-fitted with the 
best-fitting two blackbody plus power-law model 
for the earlier data, giving a x^=0.95 for 156 de- 
grees of freedom (ie having fixed all parameters ex- 
cept the normalizations). This suggests that there 
was no significant spectral variability either be- 



tween the two ASCA observations. 

2.6. Beppo-SAX observations 

RE J1034-I-396 was observed on 1997 April 18 

and 19 with the Bcppo-SAX Narrow Field Instru- 
ments, the Low Energy Concentrator Spectrome- 
ter [LEGS; 0.1-10 keV (Parmar 1997)], the Medium 
Energy Concentrator Spectrometer [MEGS; 1.3- 
10 keV (Boella et al. 1997)], the High Pressure Gas 
Scintillation Proportional Counter [HPGSPC; 4- 
60 k(N (Manzo et al. 1997)] and the Phoswich De- 
tector System [PDS; 13-200 keV (Frontera et al. 
1997)]. The source was not detected in the high 
energy instruments thus we report here the anal- 
ysis of the LEGS and MEGS only. 

The observations were performed with all three 
MEGS units and data from these wore com- 
bined together after gain equalization. The LEGS 
is operated during dark time only, therefore 
LEGS exposure times arc usually smaller than 
MEGS ones by a factor 1.5-3 (a factor of 2 for 
RE J1034 -1-396). Standard data reduction was 
performed using the SAXDAS software package 
(www.sdc.asi.it/software/saxdas). In particular, 
data are linearized and cleaned from Earth oc- 
cultation periods and unwanted periods of high 
particle background (satellite passages through 
the South Atlantic Anomaly). The LEGS, MEGS 
and PDS backgrounds arc relatively small and sta- 
ble (variations of at most 30% during the orbit) 
due to the satellite's low inclination orbit (3.95°). 
Therefore data quality depends little on screen- 
ing criteria such as Earth elevation angle, Bright 
Earth angle and magnetic cut-off rigidity (we ac- 
cumulated data for Earth elevation angles > 5° 
and magnetic cut-off rigidity > 6°). 

We extracted the LEGS spectrum from 8 arc 
min radius regions. This radius maximizes the 
signal-to-noise ratio below 1 keV in the LEGS. The 
MEGS image is not consistent with the instrument 
PSF and the spatial analysis will be presented else- 
where (Fiore et al., in preparation). We have used 
source extraction radii of 2 and 3 arc min to avoid 
as much as possible problems of confusion or the 
contribution of extended emission. The results for 
2 and 3 arc min extraction radii were consistent 
with each other and we report in the following, 
those obtained with the 3 arc min extraction ra- 
dius. 
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The LEGS and MEGS internal backgrounds de- 
pend on the position [see Ghiappetti et al. (1998), 
and the BeppoSAX Cookbook, 
www.sdc.asi.it/software/cookbook] . Background 
spectra have then been extracted from high Galac- 
tic latitude 'blank' fields from regions equal to 
the source extraction region, in detector coordi- 
nates. We have checked whether the mean level 
of the background in these observations is com- 
parable with the mean level of the background in 
the quasar observations using source free regions 
at various positions in the detectors. In the LEGS 
the mean local background is consistent with the 
'blank' field mean background. In the MEGS we 
found that the local background is higher than the 
blank field background by 7%. Wc have verified 
that the excess is constant in energy. We then 
scaled by the same amount the blank field back- 
groimd spectra before subtraction. 

Table 3 gives the observation dates, LEGS and 
MEGS exposure times and the detected count rate 
in each instrument. 

Spectral fits were performed using the XSPEG 
9.0 software package and public response matri- 
ces as from the 1997 August 31** issue. Spectra 
were always re-binned following two criteria: a) 
to sample the energy resolution of the detectors 
with four channels at all energies where possible, 
and b) to obtain at least 20 counts per energy 
channel. This allows the use of the statistic in 
determining the best fit parameters, since the dis- 
tribution in each channel can be considered Gaus- 
sian while avoiding possible systematic errors due 
to the linearity in the spacing of the original sam- 
pling. Gonstant factors have been introduced in 
the fitting models in order to take into account the 
inter-calibration systcmatics between instruments 
(Fiore et al. 1998). The expected factor between 
LEGS and MEGS is about 0.9. In the fits we then 
assume the MEGS as reference instruments and 
constrained the LEGS parameter to vary in the 
small range 0.7-1.1. The energy range used for 
the fits are: 0.1-4 keV for the LEGS (channels 11- 
400), 1.65-10 keV for the MEGS (channels 37-220) 

2.6.1. Spectral analysis 

LEGS and MEGS spectra were fitted with 
a number of models: single power-law, broken 
power-law, 2 power-laws and power-law plus disc 
black body. In all cases the emission spectrum was 



reduced at low energy by a cut-off due to neutral 
gas along the line of sight. The column density 
of the gas was always constrained to be greater 
than or equal to the Galactic value along the line 
of sight of 1.5 X 10^° cm~^. 

Table 4 gives the best fit parameters along with 
the obtained in each case. To ease the compar- 
ison of the Beppo-SAX observation with previous 
ROSAT and ASGA observations we have also fit- 
ted the single LEGS and MEGS spectra with a 
power-law model limited to the ranges 0.1-2 keV 
and 2-10 keV. 

The 90% Umit on a 6.4 keV(6.7 keV) iron Ka 
line in the BeppoSAX data is 165 eV (210 eV). 

2.6.2. Spectral curvature 

Figure 3 shows the fit with a single power- 
law: positive residuals are evident above 3-4 keV 
. The spectrum is clearly flattening above this en- 
ergy. The significance of the flattening is demon- 
strated in Figure 4 which shows the contour plot 
for the soft and hard slopes of the broken power- 
law model. The figure shows that the soft index, 
Flow is softer than the hard, Fhigh at greater than 
99% significance. 

The fits with two power-laws, a broken power- 
law and a blackbody plus a power-law all give 
a which is significantly smaller than a single 
power-law (Table 4). The blackbody plus power- 
law model produces an acceptable and a best 
fit absorbing column density lower than Galactic, 
unlike the other two models. Therefore either we 
are measiiring a small cohimn in the host galaxy 
of RE J1034-I-396, or the emission spectrum has 
a curvature also at low energy. While the former 
case is certainly possible, the possibility of a low 
energy turn down of the emission spectrum is in- 
triguing, because it may suggest a small black hole 
mass, if the low energy emission is interpreted in 
terms of optically thick disk emission. 

Indeed the soft X-ray spectrum is best fit- 
ted with two blackbodies which mimic a multi- 
temperature disc model. This result is consistent 
with the analysis of the ROSAT PSPG data by 
Puchnarewicz et al. (1995), who foimd a signifi- 
cant hardening of the 0.1-2 keV spectrum at low 
energies, when the column density was fixed at the 
Galactic value. 

The spectrum was also fitted with a bremsstrahlung 



6 



plus power-law model which gave a ^ of 75.5 
(for 94 degrees of freedom). Thus although a 
bremsstrahlung form for the soft X-ray component 
cannot be ruled out, the two-blackbody represen- 
tation does produce the lowest x^- Furthermore, 
a bremsstrahlung form for the soft X-ray compo- 
nent is not generally popular, because the implied 
extent of the optically-thin gas would be much 
larger than suggested by the variability of most 
NLSls. 

The preferred slope of the hard X-ray power- 
law continuum depends on the form assumed for 
the ultra-soft component. For the best-fitting two 
blackbody plus power-law model, the hard power- 
law is best represented by a slope a=1.2. This 
is slightly harder than that measured by Pounds, 
Done & Osborne (1995) but only slightly softer 
than that of AGN in general (where a ~1). Our 
analysis of the ASCA data suggest a slope for the 
hard X-ray power-law, a, between 0.7 and 1.3, de- 
pending on the model assumed for the soft com- 
ponent. 

Comparing the Beppo-SAX spectrum with the 
ASCA data, we find that the flux and spectral 
form of all three observations are consistent within 
90%. 

3. Model fits to the optical to X-ray con- 
tinuum 

In spite of its historical significance as one of 
the first NLSls with a strong ultra-soft component 
to be reported, RE J1034-I-396, unlike many oth- 
ers in its class, has exhibited a remarkable degree 
of flux and spectral stability. In this paper, we 
present a comprehensive review of measurements 
of RE J1034+396 in the EUV to X-rays over 6 
years, and over this period it has shown changes 
on long time-scales of up to 60% (compared to eg. 
RE J1227+164 whose flux changed by a factor of 
70). On short time-scales, no significant changes 
have been observed. Thus, in a study such as this 
where we wish to fit non-simultaneous optical to 
X-ray data of a NLSl, RE J 1034-^396 is the prime 
candidate. Obviously simultaneous optical to X- 
ray data, such as that which will be provided by 
the Optical Monitor and EPIC on XMM Newton, 
would be the ideal, but until such data are ob- 
tained, the quasi-simultaneous observations pre- 
sented here provide the best available. 



The 'nuclear' spectrum of RE J1034+396 is il- 
lustrated in Figure 5. This is made up of the point- 
source component extracted from the optical data 
(Section 2.1), the HST spectrum (Section 2.2) and 
the Beppo-SAX spectrum (Section 2.6). 

We have assumed that the optical-UV- X-ray 
emission originates primarily in a geometrically 
thin accretion disk around a rotating super- 
massive black hole. The radial temperature of 
the standard, Keplerian disk is calculated given 
the black hole mass and the accretion rate. We 
approximate the disk emission by local blackbody 
emission, which is modified by electron scattering 
in the disk atmosphere (Czerny & Elvis 1987). 
Relativistic corrections have been applied to the 
disk emission following Laor et al. (1989). 

We fit the optical to X-ray spectrum of RE J1034-h396 
using a grid of accretion disk models parameter- 
ized by a central black hole mass (M^/i), accretion 
rate (M) and inclination angle (Siemiginowska et 
al. 1995). Because there is no signature of the BBB 
in the UV spectrum of RE 1034-1-396 we also add 
a power-law component to the disk emission in 
order to describe the optical/UV continuum. We 
normalized a power law at 2 keV to 8x10^^ erg 
with a slope fixed to a = 1.15 to connect the 
near IR and hard X-ray points (see Figure 5). The 
implications of a possible physical significance of 
this component are discussed in Section 4.2. 

The fact that the spectral turn-over is detected 
in the soft X-ray band (Section 2.6.2) can be used 
to constrain the minimum black hole mass re- 
quired for the source. Assuming that the peak 
seen in the soft X-ray spectrum is related to the 
maximum effective temperature of the optically 
thick disk, then the physical radius at which this 
emission is generated is given by (Frank, King & 
Raine 1992): 



4 3GMM SRsch, 



Rn 



R-n 



49 

-^RSch 

36 



(1) 



(2) 



where T^s is the effective temperature, Rsch = 
2GMbh /(? is the Schwarzschild radius of the black 
hole and Rmax is the radius where the maximum 
temperature is reached. 

At a given effective temperature the disk con- 
tributes the most of the emission to the frequency 
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given by Sicmiginoswka & Czerny (1989) and Cz- 
erny & Elvis (1987): 

fcTeff = 1.625hu (3) 

This allow us to link the turn-over observed in the 
spectrum with the central black hole mass. Our 
approximation assumes that the disk is optic;ally 
thick and the entire spectrum is described by local 
blackbody emission. However, at high tempera- 
tures the disk photons are re-scattered by electrons 
in the atmosphere of the disk. Some photons gain 
energy and contribute to the soft X-ray spectrum. 
This modification allows to see the turn-over at 
higher frequencies, so a higher black hole mass is 
allowed. 

For a turn-over at ~ 0.25keV and an accre- 
tion rate corresponding to the Eddington luminos- 
ity, the mass derived from the above equations is 
~ IO^Mq. This mass can be larger if most of the 
observed photons are Comptonized in the disk at- 
mosphere. The UV photons can gain energy up 
to a mean value of 3kT. This is to ~ 0.25 keV for 
photons at the spectral peak of RE J 1034+396. 
If we take this temperature into account then the 
central black hole mass is estimated to be of the 
order 6 x IO^Mq, again assuming the critical ac- 
cretion rate. 

Using the optical to X-ray data presented here, 
best-fitting parameters for the black hole mass 
and accretion rate have been derived by construct- 
ing xt grids for models compared with the data 
over a range of inclination angles. In Figure 6 
we show the confidence intervals for accretion rate 
and black hole mass at inclinations [cos(inc)] of 
0.25, 0.5, 0.75 and 1.0 [where cos(inc) = 1.0 corre- 
sponds to a face-on disk]. The face-on configura- 
tion is ruled-out: note also that, in general, the 
highest accretion rates arc not consistent assump- 
tions built into the model. The results suggest 
that RE J1034+396 contains a relatively small 
black hole mass, with the best fitting value of 
~ 5 X lO^M©. The accretion rates are relatively 
high, but consistent with the model assumption. 
The best accretion rate value is found to be equal 
to 0.4 MEdd for a relatively edge-on inclination 
angle of cos (inc) =0.25. In Figures 7 and 8, we 
show some of the fits to the miclcar spectrum at 
cos(inc)=0.25 and 0.5, around the best-fitting val- 
ues. 

One problem with modeling these data, is that 



the signal to noise ratios in the optical and UV, 
and the number of available data points, are much 
higher than that in X-rays. This results in an 
imbalance in the significance of the optical/UV 
data over the X-rays in the xl statistic, despite 
the physical significance of the form of the X-ray 
spectrum. 

To establish which parameters the high en- 
ergy data alone would prefer, we have repeated 
this fitting procedure using only the X-ray data. 
The resulting xl grids are consistent with those 
for the full optical/UV/X-ray data, with a min- 
imum at Mbh = 5 X IO^Mq, M=0.4 AlEdd and 
cos(inc)=0.25, and the face-on geometry ruled out. 

As a further check on the robustness of the re- 
sults, we have also 'matched' the resolution of the 
optical/UV data to the X-rays. This was done 
by binning the spectra so that the number of opti- 
cal/UV data points per frequency decade was sim- 
ilar to that in the X-ray spectrum. Although this 
reduces the error on the measurement of the flux 
in each bin, the overall error remains dominated 
by uncertainties in the absolute flux determina- 
tion and by the effects of any spectral variabil- 
ity. Thus, the errors on the binned optical/UV 
data were also assumed to be ±20%. Again, the 
grids for the binned optical/UV data were consis- 
tent with those for the full optical/UV/X-ray data 
as above. The x^'s were generally higher at each 
point in the grids, but the minima had values of 
~1. 

4. Discussion 

The HST-FOS data and quasi-simultaneous 
0.1-10 keV spectrum measured by Beppo-SAX pro- 
vide very tight constraints on the form of what we 
assume to be the BBB (i.e. a single component 
which would normally encompass the optical/UV 
BBB and the soft X-ray excess) in RE J1034+396. 
No sign of significant emission from the BBB is 
seen in the FOS spectrum down to ~1250 A, while 
the LEGS spectrum measures a very steep and rel- 
atively strong soft X-ray component. There is also 
evidence for a flattening towards the lowest ener- 
gies in the X-ray spectrum (below '--^0.3 kcV), con- 
firming earlier ROSAT PSPC observations. Thus 
the unusually high temperature of the BBB in 
RE J 1034+396 makes the source unique for two 
reasons: (1) the BBB is shifted out of the opti- 
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cal/UV range completely; and (2) assuming no 
cold gas absorption intrinsic to RE J1034+396, 
the high energy turnover of the BBB is observed 
in soft X-rays at ~0.25 keV. 

RE J1034+396 is important in one further re- 
spect: it shows no evidence for significant variabil- 
ity on short time-scales (~days), and only mod- 
erate variability on long time-scales (months to 
years) . This is contrary to the nature of Seyferts in 
general, where sources with very soft X-ray spectra 
tend to show large amplitude variability on short 
time-scales, while harder sources vary less [see eg. 
Fiore et al. (1998) and references therein]. 

4.1. Black hole mass and accretion rate 

This unique set of observational properties, and 
the availability of quasi-simultaneous data, has en- 
abled us to place meaningful constraints on the es- 
sential parameters of RE J1034-I-396, ie. the black 
hole mass, accretion rate and inclination of the 
AD to our line of sight. The xt grids favor almost 
Eddington accretion rates (~0.4MEdd) onto a rel- 
atively low mass black hole {M^h ^ 5 x 1O''M0), 
for a disk inclined relatively edge-on (ie. 75° from 
the axis of the disk). This is consistent with the 
hypothesis of Pounds, Done & Osborne (1995), 
who suggested that RE J1034-I-396, and NLSls 
in general, were systems analogous to Galactic 
Black Hole Candidates (GBHCs). Models of GB- 
HCs predict that, in their 'high-state', the accre- 
tion rate onto the black hole is high, resulting in 
a strong soft thermal component and a relatively 
soft, fluctuating hard power-law. 

If, in the case of NLSls, the matter was accret- 
ing at high rates onto a relatively low mass black 
hole, this would also result in a relatively high 
temperature for the inner edge of the AD. Thus 
a very hot, soft X-ray component would be ob- 
served in NLSls when compared to other AGN. In 
addition, the low-mass black hole might result in 
a low- velocity broad line region (BLR), either be- 
cause the clouds were moving more slowly due to 
the weaker gravitational potential, or because the 
emitting regions were located further out (which 
in turn might be due to the weaker gravitational 
field and/or the unusual form of the incident ion- 
izing continuum). 

Our results are consistent with this idea of high 
accretion onto a relatively low mass black hole. 



There is one possible problem that should be ad- 
dressed however, and that is the inclination of 
the AD. The fits prefer an almost edge-on disk 
(~ 60 — 70°), which is not unexpected as such, 
since the probability of observing an AD in soft 
X-rays, assuming isotropic emission of the soft X- 
ray flux, is proportional to the sine of the incli- 
nation angle. This in itself is not a problem, but 
one might naively expect any molecular torus to 
lie roughly co-planar with the AD. Given that the 
opening angle of such tori is thought to be about 
45° and the optical depths implied are high [~ 
10^^ cm~^; see eg. Antonucci (1993)], this would 
effectively block out any soft X-ray or UV emis- 
sion from the disk, yet significant amounts of both 
are observed. Therefore if our modeling is correct, 
either the molecular torus is geometrically-thin, 
optically-thin, patchy or absent altogether, or the 
planes of the AD and the torus are more orthogo- 
nal than coincident. 

Boiler et al. (1997) suggested that NLSls in 
general may contain disks viewed relatively face- 
on which would contradict the results of our 
modeling. However, they based this hypothesis 
on the high levels of variability seen in many 
NLSls (and in IRAS 13224-3809 in particular) 
and RE J1034-I-396 has shown only relatively 
small changes over long time-scales. Thus while 
the variability may be related to the inclination 
of the disk, our model fits suggest that the ultra- 
soft X-ray emission, the low- velocity BLR and the 
direct or indirect link between the two, may not 
be. 

4.2. The optical/UV 'power- law' contin- 
uum 

When the FOS data are combined with the 
spectrum of the nuclear optical continuum (see 
Section 3; Fig. 5), we find a very flat, power- 
law-like component with a spectral index, a ~1, 
covermg the 8000 A to 1200 A range. There is 
no sign of any rise towards the blue in this re- 
gion (which is characteristic of the BBB), and the 
emission from the host galaxy of the AGN has 
been removed. The residual component appears 
to be non-thermal in nature. It cannot be due 
to the sum of thermal dust components because 
dust grains sublime at 10^ K and their emission 
is negligible at wavelengths short-ward of Ifim. 
Therefore, this represents the first unambiguous 
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detection of a non-thermal optical/UV continuum, 
component in any AGN. 

When modehng the optical to X-ray data, it 
was necessary to assume an underlying optical 
to X-ray power-law, in addition to the accretion 
disc component, to fit the optical/UV continuum 
and the hardening of the X-ray spectrum towards 
harder energies. The slope of this power-law was 
a=1.15, slightly softer than that in the optical/UV 
alone. The resulting deficit in the UV between 
the data and the power-law was well-modeled by 
the accretion disc component (see Figs 7 and 8). 
Thus the data and the modeling presented here, 
are consistent with the presence of a single power- 
law component which underlies the optical to X- 
ray spectrum of RE J1034-I-396. 

The existence of a single power-law component 
connecting the IR to the X-rays has been sug- 
gested before (Malkan 1984; Elvis et al. 1986; Car- 
leton et al. 1987), although the 'bare' power-law, 
such as we have observed in RE J1034-I-396, has 
not (to our knowledge) been observed, due to the 
dominance of the 'thermal' components in the op- 
tical and IR, respectively. The IR to UV contin- 
uum is generally thought to be the sum of thermal 
components from dust and the accretion disk. The 
hard X-ray power-law is considered to be physi- 
cally unrelated and produced by inverse Compton 
scattering of blackbody photons from the disk in a 
surrounding hot corona. Previous claims of a con- 
tinuous IR-to-X-ray power-law component have 
been based on correlations between IR and X-ray 
luminosities (Malkan 1984; Carleton et al. 1987) 
and extrapolations of the medium X-ray spectrum 
into the IR (Elvis ct al. 1986). In addition, Fcrnan- 
des, Sodrc & Vicira da Silva (2000) foimd evidence 
for a stable, underlying rod continuum component 
in the optical spectra of PG Quasars. The direct 
observation of the optical/UV power-law compo- 
nent in RE J1034-I-396 makes this a prime target 
for the study of broadband, non-thermal emission 
in Seyfert galaxies. 

Using ROSAT-PSPC, optical, lUE and archival 
IRAS data, Puchnarewicz ct al. (1995) also dis- 
cussed the possibility of an underlying IR to X-ray 
power-law component in RE J1034-I-396. These 
data suggested a softer slope to the power-law 
however (a ~1.4), due to the lack of any con- 
straint in hard X-rays {E >2keV) and probable 
contamination of the optical spectrum by the host 



galaxy. 

The similarity of this component to that ob- 
served in BL Lac objects led Puchnarewicz et al. 
(1995) to hypothesize that there may be some kind 
of 'mini-BL Lac' activity in RE J1034-f 396. Al- 
though subsequent optical spectropolarimetry of 
RE J 1034-1-396 revealed an upper limit on the po- 
larization of the continuum of 0.4% (Breeveld & 
Puchnarewicz 1998), the 'duty cycle' of X-ray se- 
lected BL Lac objects (i.e. the fraction of time 
spent in periods of significant variability and po- 
larization) is ~40% or less (Januzzi, Smith & El- 
ston 1994; Heidt & Wagner 1998). Thus it is pos- 
sible that RE J1034-I-396 does have a BL Lac com- 
ponent which dominates in the optical/UV, but is 
currently in a quiescent state. 

5. Conclusions 

We have presented a comprehensive study of 
the narrow-hne Seyfert 1 galaxy RE J1034+396, 
spanning the optical to X-rays and six years of 
observations. Focusing on quasi-simultaneous 
Beppo-SAX, HST, EUVE and optical data taken 
in early 1997, we have (1) measured the soft and 
hard X-ray components simultaneously for the 
first time using Beppo-SAX; (2) placed tight con- 
straints on the form of the big blue bump using 
HST, EUVE and Beppo-SAX; (3) separated out 
the nuclear component from the host galaxy in the 
optical; and (4) fitted the optical to X-ray nuclear 
continuum with the combination of an accretion 
disk and an underlying power-law, to constrain the 
black hole mass, accretion rate and inclination of 
the accretion disk in RE J1034-I-396. 

Coverage of the widest possible range in X-rays, 
particularly at soft energies, is of prime impor- 
tance for NLSls, because of their strong ultra-soft 
components. These Beppo-SAX data, which cover 
the 0.1 to 10 keV range, provide significant con- 
straints on the shape, and thus the physical na- 
ture, of the X-ray spectrum in RE J1034-I-396. 
Our analysis of the Beppo-SAX spectrum re- 
veals a significant hardening of the X-ray spec- 
trum above ~3 keV and that two blackbodies, 
with A:Toff=60 eV and 160 cV best represent the 
ultra-soft X-ray component, consistent with multi- 
temperature emission from an optically-thick ac- 
cretion disc. We do not find any compelling evi- 
dence to suggest that the hard X-ray power-law in 
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RE J 1034+396 is significantly softer than that of 
non-NLSl AGN. There is the suggestion of a low- 
energy flattening in the Beppo-SAX data, which 
was also observed in the earlier ROSAT PSPC 
spectrum. This in turn implies a high tempera- 
ture for the inner edge of the accretion disk, and 
thus a low black hole mass (which is also borne 
out by the modeling). 

The results of fitting the optical to X-ray nu- 
clear continuum have shown that the data pre- 
fer relatively high mass accretion rates (~0.3 to 
0.7 MEdd) onto a low mass black hole {Mth ~2- 
10x10^ Mq) at high inclination angles (ie. prefer- 
ably edge-on, ~ 60 — 70° away from the axis of 
the disc). The very low intrinsic columns implied 
by the X-ray fits suggest that any molecular torus 
must either be optically thin, geometrically thin, 
or lie out of the plane of the accretion disk so that 
it does not obscure our line of sight. 

The optical to X-ray nuclear continuum re- 
quires the presence of a flat, power-law-like com- 
ponent in the optical/UV. This is consistent with 
an extrapolation to the hard X-ray spectrum, al- 
though an independent origin for the hard X-ray 
power-law cannot be ruled out. The existence of 
such a component has been proposed before, but 
we believe that this is the first direct observation, 
by virtue of the unusually hot big blue bump in 
RE J1G34-I-396. It appears to be stable and unpo- 
larized and we speculate that it may be a BL Lac- 
type component which is currently in the quiescent 
stage of its duty-cycle. 
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Fig. 1. — The separate nuclear (plotted as a solid 
line) and galactic (dashed line) components of the 
optical (WHT) spectrum of RE J1034+396. 




Fig. 2.— The EUVE light-curve of 

RE J1034+396. 
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Fig. 3.— The 0.1-10 keV BeppoSAX LEGS, 

MEGS spectra of RE J1034+396 fitted with a 
power-law model reduced at low energy by Galac- 
tic absorption. Lower panel shows the ratio be- 
tween data and model. Error bars are la (68%). 
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Fig. 4. — The Tio^ — ^high contour plot for a bro- 
ken power-law fit to the LEGS and MEGS 0.1- 
10 keV data of RE J1034-I-396. Contours are 
drawn at 68%, 90% and 99% significance for two 
interesting parameters. 
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Fig. 6. — xt contour plots for the AD plus power- 
law models when compared with the optical, UV 
and X-ray data for RE J1034-I-396. 
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Fig. 5. The optical to X-ray spectrum of 
RE J1034+396, combining the WHT, JfST-FOS 
and Beppo-SAX data. The underlying optical to 
X-ray power-law continuum used in the fits is also 
shown as a dashed line. Errors are la (68%). 



Fig. 7. — AD plus power-law fits to the op- 
tical to X-ray data at an angle of inclination, 
cos(inc)=0.25. 
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Fig. 8.— AD plus power-law fits to the optical to 
X-ray data at an angle of inclination, cos(inc)=0.5. 
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Table 1 

Log of EUV and X-ray observations 



Date Instrument Observation No. Exposure time (l<s) Count rate 10 ^ count s ^ Reference 



1991 Nov 18 


R.OSAT PSPC 


RP700551 


330 


3 


(1) 


1992-1993 


EIJVE All-Sky Survey 




1.6 


2.3±0.7 


(2) 


1993 Dec 27 


EUVE DS/S 




5 


2.9±0.3 


(3) 


1994 Nov 20 


ROSAT HRI 


RH701982 


1.6 


49±2 


(4) 


1994 Nov 19 


ASCA CIS 


72020000 


26.1 


3.0±0.1 (GIS2) 


(4) 


1994 Nov 19 






25.7 


3.4±0.2 (GIS3) 


(4) 


1994 Nov 19 


ASCA SIS 


72020000 


28.2 


9.1±0.3 (SISO) 


(4) 


1994 Nov 19 






28.2 


7.6±0.2 (SISl) 


(4) 


1995 May 18 


ASCA CIS 


72020010 


10.5 


2.5±0.2 (GIS2) 


(4) 








10.3 


3.3±0.2 (GIS3) 


(4) 


1995 May 18 


ASCA SIS 


72020010 


8.9 


8.4±0.4 (SISO) 


(4) 








9.0 


6.5±0.3 (SISl) 


(4) 


1996 Nov 12 


ROSAT HRI 


RP703893 


30 


61±2 


(4) 


1997 Apr 17-18 


SAX LEGS 


50172004 


21.7 


6.9±0.2 


(4) 


1997 Apr 17-18 


SAX MEGS 


50172004 


43.2 


2.16±0.08 


(4) 


1997 Apr 14-20 


EUVE DS/S 




97.1 


3.8±0.1 


(4) 



'^(1) Puchnarcwicz et al (1995); (2) Fruscione et al. (1996); (3) Fruscione, private communication; (4) this paper; ASCA 
GIS count rates are over 1.0-10 keV, SIS count rates 0.7-10 keV 
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Table 2 
Fits to ASCA spectra 



Model 



intrinsic nn Ois or kT (eV) second BB kT (eV) 



xVdof 



72020000 Observation (1994 November) 

One power-law 0.0<0.51 

Two power-law 0.0<28.6 

Broken power-law 0.0<7.4 

Brems-|-power-law 0.0<16.7 

BB -I- power-law 0.0<3.7 

2BB -I- power-law 0.0<40 



2.04±0.09 
1.05 xlO"-'' 
o 4, +2. 08 

7.32x10-* 



4-0.45 
0.23 



2.61 



1.07x10" 



1.68x10" 
127+;^ 



4.68x10" 



6.43x10" 



5+157 
-183 

;.77xio-' 



323; 



79+0 *8 

^' ' ^ — 0.62 

3.07x10"* 



1.27"* 



74-0.20 
-0.28 

break 1.51 keV 



1.18±0.25 
5.54x10-* 



a-F0.i9 
-0.20 
6.25x10"* 



1.29+ 



r; + 0.65 
-1.25 

2.57x10"* 



0.651 



519/410 
363/408 
375/408 
367/408 
371/408 
362/406 



72020010 Observation (1995 May) 



Model 

One power-law 

Two power-law 
Broken power-law 

Brems-t-power-law 

BB -t- power-law 

2BB -t- power-law 



intrinsic nn 
0.0<2.2 



as or kT (eV) 
normalization 



an 

normalization 



37: 



56^ 



normalization 

1.98±0.13 
1.04x10"^ 

Too poorly constrained - no improvement over one power-law 

4 3, +3.62 , 73 + 0.62 

break 1.40 koV 



, +3.62 
-1.83 

2.54x10"^ 



151 



+ 117 
48 

1.16 



, + 0.60 
-0.49 

1.15x10"^ 



1.82"f 



46<99 



104+ 



1.79+ 



a+0.55 
-0.44 
1.11x10"^ 



6.05x10"* 

Too poorly constrained - no improvement over one blackbody 



151/156 

120/154 
123/154 
124/154 



'^Intrinsic nn is the column density local to the AGN, in addition to a Galactic absorbing column density 
fixed at 1.5x10^^ cm~^, and corrected for redshift. Limits (given in brackets) are 90%. Power-law indices, a, 
are defined such that Fjy oc i^"*^. Bremsstrahlung and blackbody temperatures are given in eV and have been 
corrected for redshift. Power-law and broken power-law normalizations are in units of photons keV""'" cm"'^ s~ 
at 1 keV; bremsstrahlung normalization — 3.02 X 10~^'' / {AttD'^) f neriidV, where Tie is the electron density 
in cm~'^, nj is the ion density in cm"'^ and D is the distance to the source in cm; blackbody normalization 
= Lzq/D^q, where L39 is the source luminosity in units of 10^^ erg s~^ and Dio is the distance to the source 
in units of 10 kpc. 



Table 3 

Log of Beppo-SAX observations 



Dates 


Exposure 




Rate (10-2 c/s) 






LEGS 


MEGS 


LEGS 


MEGS 


18-19 Apr 1997 


21.7ks 


43.2ks 


6.9±0.2 


2.16±0.08 
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Table 4 
Fits to Beppo-SAX data 



Model N% 



MECS 2-10 keV 

power-law 1-5(F) 



LEGS 0.1-1.5 keV 

power-law 2.3±0.5 



as or kT^ El^^^ or kT^ 



1.42±0.26 
2.45±0.20 



120.2/92 
1.13tg-^| 70.0/90 
1.08±0'.32 83.2/91 
1.50±0.26 95.1/90 
1.84±0.25 104.6/91 
1.19±0.24 65.7/89 



31.5/45 



29.9/41 



LEGS plus MEGS 0.1-10 keV 

power-law 1.75±0.35 2.27±0.14 
2.4±0.5 
1.5(F) 



broken p.l. 
broken p.l. 
bb+p.l. 
bb+p.l. 
bb+bb+p.l. 



0.5 



+0.5 



0.35 

1.5(F) 
1.5(F) 



2.54±0.16 

2.24±0.07 
kT=100±10 
kT=82t^^ 
kTi=55±10 



Eb = 2.3±0.5 
Eb = 2.5±0.6 



kT2=155t^^ 



Absorbing column density in units of 10^" cm ^; Temperatures in eV. Break 
energies in keV; (F) frozen parameter. Errors are 90%. 
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